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ABSTRACT 


The  ionospheric  contributions  to  the  non-relativistic  Doppler 

shift  at  VHF  from  satellites  above  the  ionosphere  are  considered  to  0(l/f^) 

F  zv'Z  ~  ^ 

whereof-  is  the  satellite  transmitter  frequency  and  the  vacuum  Doppler  shift 
s  f'f?  «■ 

is  considered  of .-^Ct  is  shown  that  to  0(l/f^),  the  phase  of  the 


electromagnetic  radiation  field  froSrthe-^s_atellite  cannot  be  approximated 


by  the  usual  Fermat  integral  of  geometric  optics. 


aside rat ion 


of  the  characteristics  of  an  ideal  doprpler  tracking  receiver,  boundary 
conditions  are  imposed  on  the  solution  to  the  wave  equation  containing  the 
ionosphere  electron  contribution  to  the  refractive  index  such  that  the 
solution  to  o(i/f^)  in  the  phase  is  obtained.  It  is  shown  that  to  this 
order  Fermat’s  principle  can  still  be  applied  if  the  index  of  refraction  is 
modified  by  terms  containing  gradients  of  the  electron  density.  This 
generalized  Fermat’s  principle  is  used *to  obtain  the  ionospheric  contributions 
to  the  Doppler  shift  to  0(l/f^)  at  VHF.  Upper  bounds  for  the  various  terms 
are  estimated.  It  is  shown  that:  l)  additional  terms  not  given  by  the 
geometrical  optics  approximation  are  negligible  so  that  geometrical  optics 
theory  is  valid  at  VHF  to  0(l/f^),  and  2)  except  for  very  disturbed  conditions 
in  the  ionosphere ,  the  use  of  the  two  frequency  doppler  data  to  eliminate 
the  first  order  refraction  contributi  n  should  yield  negligible  higher  order 
refraction  errors  so  long  as  the  lower  of  the  two  frequencies  is  above 
100  mc/s .  ^ 


The  principal  conclusions  of  this  paper  are  reported  in 


-  ii  - 


reference  10. 
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IONOSPHERIC  CONTRIBUTIONS  TO  THE  DOPPLER  SHIFT  AT  VHF 
FROM  NEAR  EARTH  SATELLITES 


INTRODUCTION 

Considerable  effort  is  being  spent  on  developing  methods  for  using 

accurate  measurements  of  the  radio  Doppler  shift  for  satellite  tracking  and 
/ 1  o 

geodesy. 3DJ  The  satellite  frequencies,  f  ,  currently  being  employed 
are  predominantly  in  the  50  to  500  mc/s  region  where  the  ionosphere 
contributes  significant  additions  to  the  vacuum  Doppler  shift.  Accurate 
doppler  measurements  have  utilized  the  dispersive  nature  of  the  ionosphere 
to  eliminate  the  first  order  ionospheric  refraction  contribution  which  is 
of  the  order  of  l/f  .  This  is  done  by  receiving  two  coherent  frequencies 
transmitted  from  the  satellite  (usually  a  factor  of  two  to  six  apart)  and 
combining  them  experimentally  to  eliminate  the  first  order  ionospheric 

m 

contribution  on  a  data-point-by-data-point  basis.  With  the  use  of  only 

two  frequencies  the  re.ulting  "vacuum"  Doppler  shift  still  contains  all 

contributions  that  depend  upon  higher  powers  of  (l/f  ).  In  particular, 

s 

the  Faraday  rotation  contribution  remains  as  well  as  contributions  resulting 
from  sharp  gradients  in  the  ionosphere  electron  density. 

To  this  date,  there  is  some  experimental  evidence  that  ionospheric 

contributions  higher  than  0(l/f  )  are  significant  for  frequencies  above  100 

,  (9) 

mc/s.  '  However,  because  of  the  low  powe  ’  levels  transmitted  from  those 
satellites  containing  multiple,  coherent,  frequencies,  the  noise  level  in  the 
experimental  doppler  data  is  sufficiently  high  that  combining  three  or  more 
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satellite  frequencies  to  accurately  measure  the  higher  order  contributions  have 
not  been  as  definitive  as  one  would  like. 

The  principal  objective  of  this  paper  is  to  consider  theoretically  the 
refracted  Doppler  shift  to  0(l/f^)  in  order  to  determine  whether  contributions 

\ 

higher  than  0(l/f  )  should  be  significant.  The  principal  motivation  for  this  ! 

s 

study  is  to  aid  in  the  experimental  justification  that  two  coherent  satellite  i 
frequencies  are  adequate  to  yield  accurate  Doppler  shift  data  for  tracking  and  J 


geodesy. 


The  computatic.i  of  tne  ionospheric  contributions  to  o(l/f^)  at  VHP 


frequencies  is  made  relatively  easy  because  the  transmitter  frequencies  are 
far  above  any  significant  resonance  frequencies  in  the  ionosphere.  However,  it 

O  | 

is  shown  in  Section  II  that  to  0(l/f3)  the  usual  optics  approximation  (through  '{ 

using  Fermat's  principle)  is  not  sufficient  and  consequently  the  more  difficult  I 

task  of  finding  the  actual  solution  to  Maxwell’s  equations  to  0(l/f^)  in  the  f 

s 

ih* 

phase  is  required.  r- 

J 

jt 

Section  I  considers  in  some  detail  the  expression  for  the  vacuum  Doppler 

f 

shift  when  the  ionosphere  is  assumed  to  be  absent.  In  this  section  the  character f 

I 

istics  of  an  idealized  doppler  tracking  receiver  are  considered  to  a  sufficient  j 

f 

extent  to  derive  special  boundary  conditions  on  the  vacuum  electromagnetic  | 

t\ 

field  created  by  the  satellite  transmitter.  The  resulting  analysis  yields  1 

an  expression  for  the  vacuum  field  at  an  orbitrary  point  in  space  which  can 

I 

itself  be  used  as  a  boundary  condition  for  the  refracted  electromagnetic  field. 

In  Section  II,  the  solution  to  Maxwell's  equations  for  an  arbitrary  ionospheric  ; 
electron  density  is  obtained  to  0(l/fD)  in  the  phase  for  frequencies  near  the 


*The  vacuum  Doppler  shift  is  0(f  )  so  that  contributions  of  0(l/f^)  are  of  4th 


order  in  l/f  . 

1  s 
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satellite  transmitter  frequency.  The  Fourier  synthesis  of  the  refracted  frequency 
components  is  matched  to  the  vacuum  field  expression  in  Section  III.  The  two 
fields  are  matched  in  the  region  of  the  satellite,  where  it  is  assumed  that 
the  satellite  is  sufficiently  above  the  ionosphere  that  the  electron  density 
is  negligible.  Section  III  also  presents  the  detailed  results  of  the  expansion 
of  the  phase  in  powers  of  l/fg  to  0(l/f^),  and  the  expression  for  the  refracted 
Doppler  shift  that  results.  Section  IV  presents  an  estimate  of  the  magnitude 
of  the  higher  order  terms  in  the  expression  for  the  Doppler  shift  where  it  is 
shown  that  except  possibly  for  Faraday  rotation  and  extremely  sharp  electron 
density  gradients,  the  higher  order  terms  are  most  likely  negl_gible  at  VHF. 
Appendices  I  and  II  present  some  of  the  more  tedious  analysis. 
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I.  The  Vacuum  Doppler  Shift 

Figure  1  indicates  schematically  an  ideal  doppler  tracking  receiver.  It 

is  assumed  that  the  antenna  is  circularly  polarized  and  tracks  (in  angle)  the 

satellite  to  minimize  fading  due  to  ionospheric  Faraday  rotation.  In  Figure  1, 

the  circularly  polarized  component  of  the  satellite’s  transmitter  field  that 

matches  the  polarization  of  the  antenna  is  labled  E  (t).  The  output  of  the 

s 

receiver  has  a  phase  represented  by 

03 

exp  [-  io3jj,t  +  i  «s(t)]  , 

which  is  fed  to  the  tracking  loop.  The  signal  from  the  receiver  is  mixed  with 
the  signal  from  the  VCO  to  produce  a  nearly  constant  angular  frequency  output 
at  (aXjrp  -  tDp).  This  signal  is  passed  through  an  IF- amplifier,  Bjj,(Acd)  , 
centered  at  (ox^  -  03^),  whose  bandwidth  is  Ao>and  whose  phase  response  is 
made  as  linear  as  possible.  The  phase  of  the  signal  after  passing  through  the 
IF- amplifier  is  represented  as 

03 

exp  [i  — —  A$(t)  -  i^jj,  -  03p)t] 

in  Figure  1. 

This  signal  i~  fed  to  a  phase  detector  along  with  a  reference  signal 
whose  angular  frequency  is  -  03^)  and  the  output  of  the  phase  detector  is 
presented  by 

03 

-~-A§(t) 
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Fig.  1  IDEAL  DOPPLER  TRACKING  RECEIVER 
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This  phase  difference  is  then  fed  through  a  servo  amplifier  with  response 
T(t  -  t 1 )  producing  a  phase  error  in  Figure  1  by 

CD 

8*(t). 

Tnis  phase  error  alters  the  phase  (time  dependent  frequency)  of  the  VCO  in 

CD 

such  a  way  that  the  tracking  loop  produces  as  small  a  phase  error  — —  &§(t), 
as  practicable  In  particular,  the  servo  amplifier,  R(t  -  t1),  has  enough 
prediction  and  long  enough  time  constants  that,  for  frequency  changes  typical 
of  satellite  Doppler  shifts,  the  phase  error  remains  very  small  once  the  loop 
has  "locked  on".  The  doppler  data  is  then  taken  as  the  time  derivative  of  the 
phase  of  the  output  signal  from  the  VCO. 

For  present  purposes,  it  is  sufficient  to  neglect  noise  contributions 
to  the  phase  and  to  consider  the  tracking  loop  as  ideal,  so  that 

5§(t)  =  A#(t)  =  0. 

For  this  idealized  case,  the  phase  error,  A§(t),  can  be  represented  to  a 
sufficient  approximation  by: 


M(t)  =  $g(t  -  tp)  -  §F(t  -  tD) 


where  t^  is  the  time  delay  in  the  tracking  loop  IF  amplifier  and  is  the  order 
of  the  reciprocal  of  the  IF  bandwidth,  Aco.  Consequently,  to  a  sufficient 
approximat ion : 


to 


M(t  +  tD)  =  $F(t) 


a>_ 


$g(t)  =  0 
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(1A) 


c 


c 
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and  the  experimental  Doppler  shift  is: 

f  ,  ,  to  , 

-ir  V*)  -  -  -£r- -?-■ -Ir  «.<*>  •  <“> 

Present  day  tracking  loops  have  their  IF  amplifier  bandwidth,  Ao>, 
between  approximately  5  kc/s  and  5  mc/s.  Thus  the  time  delay  in  the  instrumenta¬ 
tion  prior  to  detection  of  the  phase  error  is  no  greater  than 

*D  -  0(-ST>  S  IQ'4  see.  (2) 

Consequently,  the  doppler  tracking  receiver  must  ‘sample*  the  field  created 
by  the  satellite  transmitter  for  times  of  the  order  of  100  microseconds  around 
the  current  time,  t  .  In  particular,  components  of  the  field  radiated  from  the 
satellite  during  a  time  interval  of  about  2t^  centered  around  the  time 

R  (t  ) 

*o  "  *  ■  tD  '  v(t  )  ’ 

o 

where  Rg(tQ)  is  the  slant  range  from  station  to  satellite  and  v(t  )  is  the 
average  speed  of  propagation  from  satellite  to  tracking  station;  and  frequency 
components  within  about  l/t^  of  the  satellite  frequency  are  the  only  contributions 
to  the  transmitter  field  that  need  to  be  considered.  Since  the  satellite  moves 
a  negligible  distance  in  time  intervals  of  the  order  of  t^,  and  no  significant 
ionosphere  resonance  occur  near  the  satellite  transmitter  frequency,  the 
calculation  of  the  transverse  radiation  field  is  considerably  simplified. 

Consider  now  the  field  created  by  the  satellite  transmitter  in  the 
absence  of  the  ionosphere  and  in  a  coordinate  system  at  rest  with  respect  to 
the  tracking  station.  From  the  preceding  discussion,  only  the  appropriate 
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circular ly  polarized  component  of  the  trans/erse  radiation  field  need  be 


considered.  For  the  purposes  of  this  paper,  the  non- relativistic  expression 
* 

suffices,  and 


where 

Ey(r,t)  =  circularly  polarized  component  of  vacuum  radiation  field  from 
.satellite  transmitter  at  the  space  point,  r,  and  at  the  current 
time,  t,  in  a  coordinate  system  at  rest  with  respect  to  the 
tracking  station, 

r  (t  )  =  position  of  the  satellite  at  the  time,  t  , 

SO  o 

R(r,t  )=|?-r(t)|=  slant  range  from  the  satellite  to  the  space  point  ,r, 
o  so 

and  where 


is  negligible. 


*The  relativistic  field  can  be  obtained  by  assuming  a  spherical  outgoing  wave  in 
a  reference  system  at  rest  with  respect  to  the  satellite,  and  then  transforming 
to  the  coordinate  system  in  which  the  station  is  at  rest.  In  this  transformation, 
it  is  sufficient  to  use  as  the  relative  velocity  of  the  two  coordinate  systems: 


2  2 

The  additional  terms  of  0(v  /c  )  do  not  contribute  significantly  to  the  ionospheric 
contributions.  Consequently  the  additional  complication  of  including  relativistic 
corrections  has  been  avoided. 
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To  evaluate  Ey(cD,tQ)  ,  consider  the  space  point,  r,  to  he  placed  at 
the  antenna  of  the  tracking  station,  r^,  .  From  the  preceding  discussion,  it 
is  sufficient  to  choose  Ey((D,tQ)  such  that  the  field  E(r^,,  t  -  t^)  is  very- 
small  for  times  far  outside  the  region 


t  + 
o 


H(t0) 


-2tD^t-tD^to  + 


R(t0) 


and  for  angular  frequencies  far  away  from  to  =  2jtf  .  For  convenience, 

s  s 

E„(co,t  )  is  now  chosen  to  he 
V  o 


Ev(a,t0)  =  e3q>  *■" 


(a)  -  o>  )  +  icot  ] 

s  o 


which  yields 
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where 


r^  =  position  of  antenna  of  tracking  station, 

Rs(tQ)  =  slant  range  from  tracking  station  to  satellite. 

The  phase  of  the  tracking  filter  is  obtained  by  taking  the  expression  for  the 
phase  in  equation  (5B)  and  substituting  equation  (3)  with  v  =  c  into  the 
argument  of  Rg(to). 


T~  V*  -  V 


Consequently,  xron  equations  (l),  the  non-relativist ic  expression  for  the 
vacuum  Doppler  shift  becomes 


4fv(t)  = 


(5C) 


The  results  of  this  section,  while  not  new  in  themselves,  will  be 
used  in  the  following  sections  where  the  ionospheric  contributions  to  the  Doppler 
shift  are  evaluated  to  O(l/or)  .  In  particular,  the  same  idealized  doppler 
tracking  receiver  will  be  assumed  and  the  equation  for  the  vacuum  field, 
equation  (4),  together  with  equation  (5A)  for  E^(oj,to)  at  an  arbitrary  space 
point,  will  be  used  as  a  boundary  condition  for  the  refracted  field  when  in 
the  region  of  the  satellite. 
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II.  Refracted  Field  for  Single  Frequency  Compone  its  at  VHF 


In  this  section  the  solution  to  Maxwell’s  equations  is  presented  to 

O 

O(l/or)  in  the  phase  for  a  single  frequency  component,  cd,  near  u>  ,  which  obeys 

s 

the  radiation  condition.  The  principal  result  of  this  section  is  that  the  phase 

D 

to  G(l/ur)  cannot  he  assumed  to  he  given  hy  Fermat’s  integral 

r 


c 


S(r) 


dr*  *  §(r*)  n(r* ,to|m) 


r 

s 


where  the  path  is  taken  such  that  S(r)  is  an  extremum,  hut  can  he  represented  hy 


CD 


*  (r) 
c  s 


CD 

C 


dr’  •  £(r’)  f(r’,to|cD) 


r  (t  ) 
s'  o 


where  f(r,t_JcD)  is  given  hy  equation  (10)  of  this  section  and  the  path  is  taken 
such  that  5s(r)  is  an  extremum.  The  appendices  contribute  to  the  proof  of  this. 

Let 

p(r,t)  =  ionosphere  electron  density  where  it  is  assumed  that  it  changes 
sufficiently  slowly  with  time  that  its  change  during  the 
transmission  time  from  satellite  to  station  is  negligible, 

IL,(r)  =  Earth’s  magnetic  field  which  is  assumed  to  have  negligible  time 
dependence, 

e  =  charge  of  the  electron  in  M.K.S.  units, 
m  =  Mass  of  the  electron  in  M.K.S.  units, 

e  ,|X  =  permittivity  and  permeability  respectively  of  free  space  in 
°  0  M.K.S.  units, 
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Es(r,tojoj)  -  circularly  polarized  component  of  transverse  field  at ''the  angniift-r 

frequency,  us,  emitted  from  the  satellite  at  nearly  the  time  t  , 

o' 

£(r)  =  unit  vector  in  the  direction  of  signal  propagation, 

n(rj'k0la3)  =  equivalent  index  of  refraction  for  the  ionosphere  at  the  frequency 
os,  time  t  ,  and  space  point  r. 


The  wave  equation,  for  each  frequency  component,  Es(r,to)o>)  ,  is 


+  *eG>*0\*)  = 0  j 


(6a) 


where 


n  (r,to|o))  =  1  -  - 


(e  W 


2  .  %  s(r)  *  Hg(r) 

us2  1  f  (*V  )  - - - 

I  v  '  m'  0) 


(6b) 


and  where  the  sign  depends  upon  which  circularly  polarized  component  is 
considered,  I«t  the  solution  for  arbitrary  r  be  represented  by: 


EB(?,t0k) 


E0(?,to|a>) 

HI*,*,) 


exp 


(7A) 


where 

R(?,to)  -  -  ?s(to)| 


-  12  - 


1, 


Tr>#  Johrtt  Hopkint  Univartity 

APrill*  PMVttC*  LAMRATBRV 

Silver  Spring,  Maryland 


In  equation  (7A)  the  functions  E  (r,t  Jo>)  and  §  (r,t  jto)  are  chosen  to  he  real  and: 


*s&t0l®)  =  I  d?'  ‘  £(*')  f(r»,tja>)  , 


(7B) 


r  (t  ) 
s'  o 


where  the  path  of  the  line  integral  is  such  that  §  is  an  extremum.  So  far, 

s 

f(r,to|a>)  is  arbitrary. 

Substituting  equations  (7)  into  equations  (6),  and  then  equating  real 
and  imaginary  terms: 


7T1  -  2&  -VE  + 

o  o  2 


“  (n2  -  •  7§  )  E  =  0  , 

s  s  o  3 


m 


27  §  •  VE  -  2  -  R  •  Vo  +  E  =  0 

so  R  sos 


(8B) 


A  S(?,to) 


R(r,t  )  =  r  -  r  (t  )  ,  R  = 

0  s  0  R(r,t  ) 


(8C) 


In  Appendix  I  it  is  shown  that 


V$(r,to|ca)  =  s(r)  f(r,to|co) 


(8D) 
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Thus,  equations  (8)  can  he  written 


in  the  form: 


* 


A 

R 


VE 


-  2 


(n2  -  f2)  Eq  =  0 


(9A) 


2f&  •  VE  -  2  fR-s+EfV*&+E£-Vf  =  0  .  (9®) 

o  R  o  o 


In  Appendix  II,  it  is  shown  that 


.  ®(r,t  )  +  1’  (r,t  |to) 

S(r)  =  - °  ...  °  ■ 

■yi  + ?’  •  ?’ 


where 


?*  •  =  0  ,  ?’  =  0(-~-) 


and 


n2  -  f2  =  0(-^) 
u> 

Since 


f  =  n  +  O(-ir-)  =  1  +  0(-~ -) 
0)  o> 
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then 


|Vf|  =  0(-^-) 

CD 


Finally,  since  EQ(r,to|cD)  should  have  a  very  slow  dependence  upon  R,  assume  that: 


Ki  =  o(-i-) 

CD 


Using  the  above  relations,  it  can  now  he  seen  that: 


£  =  $  +  I'  +  OC-ijp) 

CD 

S  •  VE  =  §  •  ^E  +  Ot-ij-)  , 

CD 

$  •  £  =  i  +  o(-^-) 

CD 

V  *  £  =  +  t?  •  ?«  +  0(-~-)  , 

CD 


s  •  Vf  =  ft  •  Vf  +  O(-^JJ-) 

CD 
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Consequently,  equation  (9B)  reduces  to: 


A  *  VE  =  -  E  [R  •  ^f  +  V  •  §»  ]  +  0(-4-) 

00  a)4 

Substituting  this  equation  into  equation  (9A)  and  neglecting  terms  of  0( \J^E  )  , 

U 

'*  % 

the  following  differential  equation  results  for  f(r,t  |cd)  . 


„2,-  ^  ^  2/-  .  1  n  ,  C2  [R  •  Vf  +  V  •  I»  ]  ,  n/  1  % 

f  (r,t  |(d)  =  n  (r,t  Ico)  +  — ^ - -jj -  +  0(~ 5“) 

c a  w 


However,  since. 


Vf  =Vn  +  0  (-V) 


£0 


the  expression  for  f(r,to|co)  finally  reduces  to: 


,2  A 


f(r,tJco)  =  n(r,tn|a>)+  — - —  +  O(-^g-)  ,  (10) 


2a) 


0) 


Prom  equations  (7)  and  (10),  it  can  be  seen  that  to  0( — jj— )  the  phase  f 

co 

of  the  field  is  given  through  integrals  of  the  Fermat  type  if  the  integrand 
contains  the  function,  f(r,t  'to),  instead  of  the  index  of  refraction,  n(r,tJco).j 
To  0(  ■■■)—•)  ,  the  difference  between  f(r,t  |co)  and  the  index  of  refraction  is 


£0 


2u) 


R  •  Vn  +  V  •  §' 
R 
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The  first  term  depends  upon  the  gradient  of  the  index  of  refraction  along  the 

geometric  path.  From  Appendix  II,  the  vector,  §’  (r,tQJaj),  is  the  gradient  of 

the  deviation  of  the  optical  path  from  the  geometric  path  taken  in  the  direction 

of  the  geometric  path.  Therefore,  the  second  term  in  the  difference  between  f 

and  n  is  also  proportional  to  the  spatial  change  of  the  index  of  refraction. 

Consequently,  the  ’corrected  index  of  refraction’  is  proportional  to  the  spatial 

change  in  the  index  of  refraction  and  can  become  large  if  the  electron  density 

* 

in  the  ionosphere  is  not  a  slowly  varying  function  of  position. 

Substituting  equation  (10)  into  equation  (7B),  evaluating  the 
’corrected  optical  path’,  and  then  performing  the  integral  along  this  path, 

— » 
r 

Lr’  •  £(r’)  n(r’,to|co)  + — ^ 

2a> 

r  (t  ) 
s'  o' 

However,  it  is  shown  in  Appendix  II  that  th 
geometrical  path  by  an  amount,  §(r,to|co)  = 

geometrical  path  can  be  used  in  integrating  the  corrections  to  the  index  of 


te*  •  £(?•)  R  '  ^’n  +  o(; 


r  (t  ) 

S '  O 


e  optical  path  deviates  from  the 
0( — | — )  .  Thus,  to  0(— — g— )  the 


*It  should  be  noted  that  equation  (10)  was  derived  assuming  the  satellite  frequen 
is  significantly  above  any  ionospheric  resonance  frequency.  This  expression  is 
not  valid  for  frequencies  near  to  or  at  any  of  the  ionosphere  resonance  frequenc 
Equation  (10)  should  not  be  used,  for  example,  in  the  consideration  of  ionosonde 
data.  See,  for  example,  reference  (6). 
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geometric  path.  Prom  Appendix  II ,  the  vector,  §*  (r,tQ|oj),  is  the  gradient  of 

the  deviation  of  the  optical  patn  from  the  geometric  path  taken  in  the  direction 

of  the  geometric  path.  Therefore,  the  second  term  in  the  difference  between  f 

and  n  is  also  proportional  to  the  spatial  change  of  the  index  of  refraction. 

Consequently,  the  ’corrected  index  of  refraction’  is  proportional  to  the  spatial 

change  in  the  index  of  refraction  and  can  become  large  if  the  electron  density 

'ft 

in  the  ionosphere  is  not  a  slowly  varying  function  of  position. 

Substituting  equation  (10)  into  equation  (7B) ,  evaluating  the 
’corrected  optical  path’,  and  then  performing  the  integral  along  this  path. 


V*n  +  V’  -  I’ 
R’ 


However,  it  is  shown  in  Appendix  II  that  the  optical  path  deviates  from  the 

geometrical  path  by  an  amount,  §(r,tja>)  =  0( — - — ).  Thus,  to  0( — —)  the 

0)  0) 

geometrical  path  can  be  used  in  integrating  the  corrections  to  the  index  of 


*It  should  be  noted  that  equation  (lO)  was  derived  assuming  the  satellite  frequency 
is  significantly  above  any  ionospheric  resonance  frequency.  This  expression  is 
not  valid  for  frequencies  near  to  or  at  any  of  the  ionosphere  resonance  frequencies 
Equation  (10)  should  not  be  used,  for  example,  in  the  consideration  of  ionosonde 
data.  See,  for  example,  reference  (6). 
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refraction  in  equation  (lG),  and  the  phase 


6A»vu  wj  • 


CO 


c 


CD 

C 


r_ 

Jdr 1  •  £(r' )n(r*  ,tQ{cD) 


r  (t  ) 
s'  o 


A 


2cd 


dr*  •  R 


r  (t  ) 
s'  o 


R  •  V’n(r*,t  |cd)  +  V*  •  §’ (r’,t  |co) 


R(r»,to) 


+  0(-“-), 

ay 


(11) 


where  the  usual  optical  path  can  now  he  used  in  the  first  integral  in  equation  (ll' 
Equations  (10)  and  (ll)  are  the  principal  results  of  this  section.  If 
the  usual  form  of  Fermat’s  principle  is  used,  the  phase  is  given  by 


CD 


s(?)  =  -f- 


dr*  *  §(r’)  n(r’,to|cD) 

r  (t  ) 
s'  o' 


y 

es. 


where  the  path,  s(r),  is  taken  such  that  this  integral  is  an  extremum.  From 
equation  (ll),  it  can  be  seen  that 
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03 


c  (4s-S^ 


203 


27 

Jar'  •  $(r*,to)  — 


V‘n  +  V' 


R(v.t  ) 
'  l  o' 


o(-V) 


o> 


r  (t  ) 
s'  o 


and  the  right  hand  side  of  this  equation  must  he  negligible  before  the  usual 
optical  approximation  for  the  ionosphere  can  be  used  if  results  to  0(— ' ^-)  are 
desired. 


03“ 
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III.  Expression  for  the  Refracted  Doppler  Shift 


In  this  section  the  single  frequency  components  are  combined  to  yield 
the  expression  for  the  field.  The  evaluation  of  the  resulting  Fourier  intet ..  al 
is  performed  in  the  following  way. 

1.  Asstime  that  the  satellite  is  sufficiently  above  the  ionosphere  that  in  a 
region  near  the  satellite  the  electron  density  is  negligible, 

2.  match  the  expression  for  the  refracted  field  to  the  vacuum  field  at  space 
points  close  to  the  satellite, 

3-  assume  that  when  the  space  point  is  placed  at  the  tracking  station  antenna, 
no  singularities  are  introduced  into  the  amplitude  of  each  frequency 
component  for  frequencies  near  the  satellite  frequency,  and  then 

4.  evaluate  the  integral  over  the  frequency  domain  by  the  method  of 
steepest  descents. 

O 

The  resulting  expression  for  the  phase  is  then  expanded  to  0(l/o>  J)  and  the  time 


derivative  taken  to  yield  the  expression  for  the  refracted  Doppler  shift  to 

0(l/cDg^)  of  the  vacuum  Doppler  shift  itself.  Appendix  II  contains  a  detailed 

calculation  of  the  extremum  path  for  the  phase  function,  §  ,  ,  which  is  required 

s 

for  the  final  evaluation  of  the  phase  in  powers  of  l/co  . 


*  *r 

From  equation  (7A) ,  the  expression  for  the  field  at  an  arbitral 


'Ovnr  cno  r\e. 


point  is: 


E  (r,t  It)  = 
s'  3  o'  ' 


E  (r;t  |cd) 

- - -  exj{-icoCt 


-  ~ ~  $  (?,t  1(d)]]  . 
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Now  assume  that  the  satellite  is  sufficiently  above  the  ionosphere  that  a  region 
in  the  vicinity  of  the  satellite  can  he  found  where  the  far  field  pattern  of  the 
radiation  field  is  still  essentially  the  vacuum  field  and  the  electron  density  is 
negligible.  For  space  points  r  ,  in  this  region  the  refraction  index  is  unity 
and  equation  (12)  reduces  to: 


CO 


R(r,to)  -ji 


E  (r,t  ), 
v  3  o  3 


where 


R(r,t  )  «  R  (t  )  =  |r  -  r  (t  )j 
'  3  o  s'  o  1  T  s'  o  1 

From  equations  (b)  and  (5A.),  it  can  be  seen  that  within  this  region, 

-  ^r-  (« -  «.>2  * 

R(r,t  )  «  R  (t  ) 

'  3  o  s'  o 

Let  the  field  point  now  be  placed  at  the  antenna  of  the  doppler  tracking 

station,  r  =  ,  and  assume  that  for  frequencies  in  the  complex  plane  near  the 

(real)  satellite  frequency,  co  ,  no  singularities  or  branch  points  appear  in 

s 

El^T,t0la)^  Where: 


E0(r,tJcD)  =  exP 
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E  (r  ,t  Ico)  =  E  (r  ,t  Ico)  exp 
o  t  o'  1  T  o' 


j 


It  is  shown  in  Appendix  II  that  no  singularities  or  branch  points  appear 

in  the  phase  for  o  near  oo  ,  and  therefore  the  integral  in  equation  (if)  for 

s 

the  radiation  field  can  be  evaluated  by  the  method  of  steepest  descents. 
Equation  (12)  then  reduces  to: 


W 


exp 


1  [t  -  t  -  —  (co  $  )] 

dco  s  s' 


“d 


o  c 


V  2"  V5  ' 


exp  jtog[t  -  tQ  -  j  • 


(13) 


The  time  interval,  t^,  is  the  order  of  10”  sec.  and  in  this  time 
the  satellite  moves  a  negligible  amount.  Thus,  in  the  phase,  the  time,  tQ, 
can  be  identified  with  the  time,  t,  by  the  expression 


*0  =  *  -  "T-gl-  C“sls^T’tK)]  •  (l4A) 

s 

The  Doppler  shift  is  proportional  to  the  time  derivative  of  the  phase,  and 
thus  the  refracted  Doppler  shift  is  represented  by: 
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Af(t)  = 


1  %  d 


2rr  c 


W  VV^K’ 


(14B) 


It  now  remains  to  evaluate  the  function  §  in  detail. 

s 

In  Appendix  II,  it  is  shown  that  the  optical  path  deviates  from  the 

2 

geometrical  path  "by  an  amount  of  0(l/o3  ).  Thus,  from  equation  (6b), 

s 


n(?,tol<os)  =1-  -i-  p(r,t0)  ±  -i-  -si-  -j£-  p(?,to)  gs(to)  •  Hg© 


2  2mc 

03  o 

s 


os  2m  o 
s 


1  e 


4  fl2  2 
03  om  e 
s  o 


P2(r,t  )  +  0( — ~-) 


03 


ft  ft  ,  rT~  *,(*„> 

S  °  l%-?s<VI 


Substituting  this  expression  into  equation  (10),  and  defining: 


a2(r.to)  = 


2me 


p(?,to)  , 


(15A) 


■36*0*  ■  *  “i2-  ’ 

zttl  O 


(15B) 
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Vr*V  =  " 


2  2  R  (t  )  •  vp(r,t  ) 
r  p.  s'  n  3  <v 


8m\2  p2(r>V+-^0  |?-rs(tQ)| 


2  2~*  - 1  /—  ,  » 
c  to  V  '  §  (r,t  cd  ) 
s  '  ’  o1  s 


the  expression  for  f(r,to|  cDg)  "becomes: 


iL  a,  (r,t  )  , 

f(?,t  l»s)  =  i  +  £  ^-1r^-  +  0(-V) 

m  CD  CD_ 


i  m  CD 

k=2  s 


(16) 


Also  in  Appendix  II,  the  integral  of  this  function  over  the  optical 

path  is  given  to  0(i/cd  )  hy: 

s 


=  d^‘  Rc(t  J  Cf(r*  ,tJcoJ  -  1?  (r  ,tJcDc)|2]  +  0(-~-)  , 


S'  o 


o  *  s /  2 


o'  s 


CD 


r  (t  ) 
s'’  o 


where  ^s(tQ)  •  §  (r,tQ|cDs)  =  0,  and  each  of  the  two  components  of  §  = 


are  given  hy: 


(17A) 
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k  =  1,2  . 


r 

f  -+*  A 

The  symbol,  j  dr 

r  (t  ) 
s'  o 


Rs(to)  denotes  taking  the  integral  along  the  geometric 


path.  The  derivatives 


da2(?  ,tQ) 

N  ^ 

OSk 


§  =  0 


k  =  1,2, 


denote  taking  the  gradient  of  a?(r,to)  in  two  orthogonal  directions,  both 
orthogonal  to  the  direction  of  the  geometric  path,  and  evaluated  on  the  geometric 
path. 
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Finally,  substituting  equation  (l6)  into  equation  (l/A)  and  letting 


y  w = 


(air)2  J 


T 

r 

I  -»  A  .  . 

dr  •  R  (t  )  a0(r  ,t  ) 
s'  Qs  2 '  J  o 


(18A) 


r  (t  ) 
s  o 


V  W  = 


(2n): 


T 


t  Ay  .  . 

dr  •  R(to)  a3(r  ,tQ) 


(l8B) 


r  (t  ) 
s'  o 


*T 


vw = 


(2TT)"  j 


r 

r  2 

(to  •  S(t  ) 

f-*1  j_  \  ®s  i  i?1  z-*’  j.  i  1 2 

a4(r  'V  -  2  1  5  (r  ^okl 

9 

— 

r  (t  ) 
s'  o7 


(180) 


it  can  be  seen  that  the  phase  is  given  by: 


J*' 


§  (r,t  jco  )  =  R  (t  )  +  \  ( - ^T— |  #  (R  ,t  )  +  0  ( — — jrH 

s'  T  ol  s7  s'  o7  ^  f  k  J  k'  s;  o7  5  J 

k=2 


(19) 


where  f  =  o  /2rr  has  been  used, 
s  s7 
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To  order  v/c  (non-relativist ic  approximation),  the  time  derivative 
of  the  phase (to  give  the  Doppler  shift) is  given  by  taking  the  time  derivative 
of  the  phase  with  respect  to  t  and  then  substituting  the  value  of  tQ  given 
by  equation  (l4A) .  letting 


- 


5 

*k+l 


(R  ,t  ) 


(20A) 


the  expression  for  the  refracte*  Doppler  shift  becomes 


df(t)  =  4fv(to)  + 


> 


(20B) 


where,  from  equations  (l4A)  and  (19) ,  to  0(v/c): 


t 

o 


=  t  - 


R  (t) 
sN  ' 


(k  -  1) 

k 


+  0/ 


/  1 


n  5 

Cl 

s 


(20c) 


Equations  (20),  together  with  equations  (15),  (17B)  and  (18),  provide  t.*e  required 
expressions  for  the  non-relativistic  refracted  Doppler  shift.  The  expressions 
are  valid  to  0(1/f  ~>)  or  to  of  the  vacuum  Doppler  shift  itself.  Section 

IV  contains  a  summary  of  these  results  and  some  numerical  up^er  bounds  to  aid  in 
estimating  the  significance  of  the  various  terms. 
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IV.  Summary  and  Conclusions 

! 

I 

The  equations  given  in  Section  III  for  the  refracted  Doppler  shift 

i 

are  written  in  a  notation  especially  chosen  to  display  those  variables  important 
in  matching  boundary  conditions.  This  previous  notation  is  somewhat  awkward 
,  for  purposes  of  discussion  and  the  final  equations  of  Section  III  will  now  be 
rewritten  in  a  slightly  different  notation.  Figure  2  schematically  chows  the 
geometry  of  the  satellite  motion  during  a  small  time  increment.  At  ,  relative 

* 

to  the  tracking  station.  For  convenience,  the  origin  of  the  new  coordinate 

system  is  taken  at  the  tracking  station  antenna.  From  Figure  2,  it  can  be  seen 

j  that  the  velocity  vector,  §  (t  )  ,  lies  in  the  plane  defined  by  the  position  and 

— * 

velocity  vectors  of  the  satellite,  rs(tQ)  and  rg(tQ)  respectively,  and  is  normal 

to  rs(tQ;  .  Finally,  it  can  be  seen  that  the  new  unit  vector,  rg  ,  is  just  the 

A  .  . 

negative  of  the  unit  vector,  Rg(to;,  used  in  the  previous  sections 

With  this  notation,  for  an  arbitrary  function,  g(r,t  )  ,  and  r  =  rr  : 

o  s 


Also  letting 

o  2m®  j 

=  a>E  )  , 

€ 

*  =  1,2, 

the  various  expression  for  the  refracted  Doppler  shift  can  be  written  in  the 


following  way. 


IX  KTHi  xonm  vtmunr 
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MX  IWO  *u,WO 


TRACKING  STATIOK: 


Fig.  2  GE0M3TRY  FOR  SMALL  CHANGE  IN  t 
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Af  (t  )  =  Af  (t  )  + 

'  O  V  o 


where 


Af,  (t  ) 
k'  o 


+  0 


k=l 


(21A) 


-fv(t  )  =  vacuum  Doppler  shift. 


(21B) 


Af  (t  ) 
k'  o 


"s  & 


^ lr 4**1  >  k  “ 


dt  k+1  '  s'  o 
o 


(21C) 


t  =  t  - 


r  (t  ) 
s'  o 


* 


(22A) 


v  =  C 
P 


1  - 


i  VvV 


r  (t  ) 
s  o 


+  0 


_2_\ 


J  ’ 


(22B) 


r  (t  ) 
s'  oy 


$0(r  ,t  )  = 
2  s 3  o 


(2tt) 


2  2  me 


dr*p(r  r* ,t  ) 
s  '  o 


(23A) 


r  (t  ) 
s'  o 


VvV 


,3  U. 


r 


(Sny  2mS  C'o 


J 


dr'p(?Er’,to)  £g 


^(V')  ,  (23B) 


*If  the  relativistic  expression  is  desired,  the  relativistic  vacuum  doppler 
expression  should  he  Af  (t0)  and  equation  (22A)  altered  to  he  the  corresponding 
relativistic  expression.  The  refraction  correction  terms,  Af^(t0)  ,  are 
sufficiently  small  that  the  non-relativistic  expressions  are  still  valid. 
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r  (t  ) 
s'  o 


$,,(r  ,t  )  =  - 
4'  S ’  o' 


8m2e  2 
o 


fdr’[p2(r  r». 


1  r'^t  ) 
s  ’  o 


-  V<?*rto>  - 


V(?’,t0)] 


(23C) 


r  (t  ) 
s'  o 


2  2 
CV 


"  r 


r — ; —  p(r  r*  ,t  )  + 
L*  s  ’o' 
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~§r“  Vv*  +  ^to)  *  “SIT-  Vrsr‘  + 
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r  (t  ) 
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p(r  r*  +  ?,t  ) 


(23D) 
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|  In  these  expressions,  the  following  identity  can  he  used. 

! 


! 


r  (t  ) 

sv  o 


dt 


f 


j 

r 


dr'  h(rr*,to) 


dr  (t  ) 
s'  o' 

dt 


h(r  (t  ),t  )  + 
'  s'  o’  o 


(2b) 


Sh 

W' 

o 


where,  hy  assumption  number  1  of  Section  III,  the  first  term  is  zero  when  the 


function  h(rs(tQ) ,tQ)  is  proportional  to  the  electron  density  or  its  gradients 
evaluated  at  rs(tQ).  Except  possibly  in  Af^t  ),  the  explicit  time  dependence 

Vi_ 

of  the  electron  density  can  be  neglected,  so  that  usually  (— ^ — )  can  be 

o 

neglected  in  equation  ( 2b ). 

Equations  (2l)  through  (23)  are  the  non-relativistic  expressions  to 

0(l/f?)  for  the  Doppler  shift  corresponding  to  a  circularly  polarized  receiving 
s 

antenna.  Admittedly  they  are  rather  complicated.  As  expected,  the  lead  term  is 

the  vacuum  Doppler  shift,  evaluated  at  the  satellite  position  corresponding  to 

the  retarded  time,  t  -  rs(t0)/v  ,  where  v^  is  the  average  signal  velocity  to 

0(1/1?).  Higher  corrections  to  the  signal  velocity  are  negligible  when 
s 

considering  the  non-relativistic  Doppler  shift. 

The  first  order  refraction  term,  Af^(to)/fg,  is  the  term  that  is 
eliminated  when  two-frequem.y  doppler  data  is  used.  From  equation  (23A)  it  can 
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be  seen  that  this  first  order  term  is  proportional  to  the  time  rate  of  change 
of  the  total  number  of  electrons  in  a  tube  of  unit  cross-section  extending 
from  the  tracking  station  to  the  satellite.  Figure  3  shows  a  typical  set  of 

experimental  data  taken  from  the  satellite  i960  ETA  1.  On  rare  occasions, 

8  2 

Af^(t  )  x  10~°  has  been  observed  to  be  as  large  as  10  (cps)  . 

2 

The  second  order  contribution  Af  (t  )/f  ,  can  be  seen  from  equation  (23B) 

o  s 

to  correspond  to  the  same  term  in  the  ionospheric  refractive  index  that  gives  rise 
to  Faraday  rotation  since  its  sign  depends  upon  whether  the  antenna  is  left  or 
right  circularly  polarized.  From  equations  (23A)  and  (23B)_,  it  can  be  seen 
that  Afg(to)/fs  has  an  upper  bound  which  is  approximately 


“2 

/  e 

■*<& : 

Afi 

„2 

~  2nm 

■r 

f 

f 

s 

~s 

’  S  i 

(25) 


Choosing  an  upper  bound  for  Af^ 


x  10 


■•8 


of  10  (cps) 


2 


Afg  x  10 


-16 


<  O.llf  (cps)' 


2 

or  at  fg  =  100  mc/s,  Af f ^  <  0.l4  cps.  Under  normal  ionosphere  conditions ,  the 
magnitude  of  this  term  should  be  a  factor  of  three  or  four  smaller. 
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Fig.  3  FIRST  ~RDEH  REFRACTION  DATA  FROM  SATELLITE  i960  ETA  1 
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From  equation  (23C)  the  third  order  contribution  is  composed  of  four 
distinct  terms. 


V(vV  =  • 


*4,2VW 


'4,3VW 
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2  2 
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s  ’  o 


(2n)4  4ffie0 


r  (t  ) 
s'  o 


J dr*  [^(r*  ,to)  +  T)22(rSto)] 


r  (t  ) 
s'  o' 


dr* 

(r  (t  )  -  r1 ) 
'  s'  o 


(2  6k) 


(26B) 


d  p(rsr’ ,tQ)  ,  (26C) 


dr’ 


s  o 
r 


a,.  :  **i  . 

:  i 


/  •  s  ’  o 

i 

o 


(26D) 


5=0  , 


where  the  1L(r,t  )  are  given  by  equation  (23D). 

X/  o 
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The  first  term  arises  from  the  fourth  order  term  in  the  expansion  of  the  refractive 
index  in  powers  of  l/fg.  The  second  term  represents  the  contribution  to  0(l/f^) 
of  the  difference  between  the  optical  path  and  the  geometric  path.  The  last  two 
terms  of  equations  (26)  represent  the  additional  terms  obtained  from  the  solution 
of  the  wave  equation  that  are  not  obtained  when  the  geometrical  optics  approxi¬ 
mation  is  used.  Crude  upper  bounds  to  the  magnitudes  of  their  contributions 
to  the  refracted  Doppler  shift  are  now  given. 

Prom  equation  (24),  neglecting  the  explicit  time  derivative. 


at 


(2  it)1 


^  frt  t  f*  .  u.  \  W 
o  J 


*  Vp(rsr‘,to) 


J 


< 


(ary 


2 

e  p  max. 
2me 


(2tt)‘ 


2me 


s 

dr* 


VP 


<  f 


I  d* 

.  2 


max  dt 
i  o 


f2  • 
max 


(*t)‘ 


2 

e  p  max 
2me  > 


where  is  the  maximum  plasma  resonance  frequency  in  the  ionosphere; 


f  „  <15  mc/s. 

max  ~  ' 
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From  equation  (21C),  using  an  upper  bound  for  Af^(tQ)  x  10"°  of  10  (cps)  , 


I  Af _  _ (t  ) I  <  f2  I  Af. (t  ) I  ,  f 
1  3>1  o'  ~  max  1  lv  oy|  >  \ 


,2 

max 


e  p  max 


(2TT  f 


(2T) 


or 


|Af3A(tQ)|  x  10‘24  <  .23  (cps)^ 


From  equations  (24)  and  (26),  it  can  be  seen  that  the  remainder  of 
the  terms  depend  upon  various  derivatives  of  the  electron  density.  To  obtain 
approximate  upper  bounds,  it  is  assumed  that  over  a  distance  of  about  one 
kilometer,  the  electron  density  can  change  by  no  more  than  its  maximum  value 
and  that  its  second  derivative  can  be  no  larger  than  ten  times  the  maximum 
electron  density. 

Accordingly,  assume: 


|vp(^r,t  )|  <  p  5(r  -  r  )  , 


j72p(^r,t  )|  <  10  p  6(r  -  r  ) 
'  ;  o' 1  ~  Kmax  N  g 


(28a) 


(28B) 
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From  equations  (24)  and  (28A) } 


r  (t  ) 
s'  o 
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1 7]  - 1  <  p  !  {dr1  6(r‘ 
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r  (t  ) 
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r  (t  ) 
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j 


<  P, 


max 


1  -  S(r  ) 
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r  -  r 

t  g 
1  “  r  (t  j 
s'  o 


l  =  1,2, 


where  S(r  )  is  the  unit  step  function  at  r  =  r  .  Bins, 
g  g 


f 


dr*  I  T^  +  Tig 


<  p  r  (t) 
~  'max  s' 


g 


s  L 


1  - 


g 


s  J 


2  rs<V 


max 


so  that 


Val  -Vs-  ^ 


Substituting  this  expression  into  (21C);  and  using  the  fact  that  the  vacuum 
Doppler  shift  at  100  me/s  is  not  greater  than  about  2500  eps, 
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f2 

max 

T~ 


i  4fv<VI 


; 


(29) 


!Af3,2(VI  X  10"24S  -OSCcps)4 

Employing  similar  analysis  using  equations  (28),  it  can  "be  shown  that 

I «3<3(t0) |  x  lCf  24  <  0[l0'5(cps)4]  , 

|Af3^(to)|  X  10"84  <  0[10J‘  (cps)4] 

Therefore ,  unless  the  satellite  is  below  VHF  (near  the  plasma  resonance 

frequencies) ,  the  two  terms  not  obtainable  from  the  geometrical  optics 

approximation  are  negligible.  The  fact  that  these  two  additional  terms  are 

negligible  affords  a  proof  that  for  each  circularly  polarized  component  of 

the  received  signal  from  the  satellite,  the  Doppler  shift  at  VHP  is  given  to 

0(l/f^)  to  sufficient  accuracy  by 
s 

Af  =  $(t)  ,  $  =  Jnds  (30A) 

where  the  integration  path  is  the  optical  path  which  makes  5  an  extremum.  From 

equations  (2l),  (22),  and  (23),  it  can  be  seen  that  to  0(l/f2), 

s 


Af  =  -fr  fndr 
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where  this  integral  is  taken  along  the  geometric  path. 

Table  I  summarizes  the  results  pertaining  to  upper  hounds  on  the 

various  contributions  to  the  refracted  Doppler  shift.  Assuming  that  these 

bounds  are  realistic,  the  second  and  third  order  refraction  contributions  may 

be  significantly  large.  It  is  believed  that  the  upper  bound  for  Af2(tQ)  is 

correct  to  within  a  factor  of  two,  so  that  if  the  ionosphere  is  unusually  dense 

and  the  propagation  path  is  parallel  to  the  earth’s  magnetic  field,  the  Faraday 

contribution  should  not  be  ignored  at  the  lower  end  of  the  VHF  region.  The 

estimated  upper  bound  for  Af  (t  )  may  be  in  error  by  a  factor  of  five  since 

o ,  x  o 

it  is  dependent  upon  the  square  of  the  electron  density. 

The  upper  bound  for  Af  (t  )  is  very  difficult  to  estimate.  It  is 

5 ){-  o 

relatively  easy  to  show  that  if  the  ionosphere  electron  density  contains  no 

sharp  gradients,  Af  (t  )  is  negligible.  Consequently,  only  unusually  large 

5)^-0 

gradients  can  contribute  significantly,  and  equations  (28)  v,ere  assumed  to  aid 
in  establishing  upper  bounds  to  the  contributions  from  regions  of  sharp  density 
gradients.  However,  equations  (28)  assume  that  only  a  single  region  contributes 
at  any  one  time.  When  the  ionosphere  is  unusually  disturbed,  for  example  during 
a  severe  magnetic  storm  or  in  auroras,  there  could  be  many  such  regions 
contributing  along  the  transmission  path.  If  these  should  constructively 
contribute  (have  the  same  sign)  the  contributions  from  density  gradients  could 
be  considerably  higher.  Such  cases  should  be  rare,  however,  since  there  is 
the  probability  that  if  there  are  many  regions  of  large  gradients  they  would 
contribute  with  random  signs  and  partially  cancel. 

Table  I  contains  a  summary  of  the  estimated  maxima.  Except  for  the 
vacuum  term,  typical  values  are  probably  a  factor  of  three  to  ten  smaller.  If 
very  accurate  measurements  of  the  vacuum  Doppler  shift  are  required,  only  time 
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periods  when  the  ionosphere  is  normal  should  "be  used.  For  such  carefully 
selected  periods,  contributions  higher  than  first  order  should  be  negligible. 

Table  II  presents  estimated  maxima  of  the  second  and  third  order 
contributions  to  the  ’vacuum’  Doppler  shift  when  dual  frequency  data  are 
combined  to  eliminate  the  first  order  refraction  contribution.  Letting  the 
two  transmitter  frequencies  be  f^  and  f^  and  the  corresponding  experimental 
data  be  Af  and  Af_  ,  the  combined  Doppler  shift  becomes: 


where  f  is  a  chosen  reference  frequency.  Table  II  lists  three  typical 
combinations  where  the  entries  in  the  table  have  been  converted  from  frequency 
shift  to  errors  in  the  range  rate  in  meters  per  second.  From  this  table,  it 
can  be  seen  that  if  the  lower  of  the  two  frequencies  is  not  less  than  lOOmc/s 
and  the  ionosphere  is  not  extremely  disturbed,  the  second  and  third  order 
refraction  contributions,  should  be  negligible. 
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periods  when  the  ionosphere  is  normal  should  he  used.  For  such  carefully 
selected  periods,  contributions  higher  than  first  order  should  be  negligible. 

Table  II  presents  estimated  maxima  of  the  second  and  third  order 
contributions  to  the  1  vacuum'  Doppler  shift  when  dual  frequency  data  are 
combined  to  eliminate  the  first  order  refraction  contribution.  Letting  the 
two  transmitter  frequencies  be  f^  and  f^  and  the  corresponding  experimental 
data  be  Af  and  Af_  ,  the  combined  Doppler  shift  becomes: 


•i 

v 

!' 


where  f  is  a  chosen  reference  frequency.  Table  II  lists  three  typical 
combinations  where  the  entries  in  the  table  have  been  converted  from  frequency 
shift  to  errors  in  the  range  rate  in  meters  per  second.  From  this  table,  it 
can  be  seen  that  if  the  lower  of  the  two  frequencies  is  not  less  than  lOOmc/s 
and  the  ionosphere  is  not  extremely  disturbed,  the  second  and  third  order 
refraction  contributions  should  be  negligible. 
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APPENDIX  I.  Proof  of  Equation  (8D). 


While  well  known,  the  author  has  not  found  in  the  current  literature 
a  detailed  proof  of  the  fact  that  if  the  line  integral, 


(Ai.l) 


is  taken  over  that  path  going  from  the  point  rQ  to  the  point  r^  which  makes 
§  an  extremum,  then 

A. 

*  *  =  f(?1)  ,  (ai.2) 

where  s(r^)  is  the  unit  vector  tangent  to  the  path  of  integration  at  the 
point  r^  .  Equation  (Al.l)  is  given  as  equation  (8d)  in  the  main  body  of 
the  paper.  A  proof  is  given  in  this  appendix  for  the  benefit  of  the  reader. 

Let  a  parameter,  a  ,  vary  from  zero  to  one  when  following  the  path 

~ ~ 

from  rQ  to  r^  ,  and  let  the  path  he  specified  by: 

r  =  r(a,rx)  =  x{o,r^)  ,  yCa,^)  ,  z(a,r1) 


*This  theorum  forms  the  basis  of  bhe  transition  from  electromagnetic  theory  to 
geometrical  optics.  See,  for  example,  references  (5),  (7) ,  and  (8). 
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Writing 


r>  =  j£- 

da  ' 


the  line  integral  can  he  written  in  the  form 


§  =  /da  L(r',r,a)  , 


L(r’,r,a)  =  [x*2  +  y’2  +  z*2]  f(x,y,z) 


(A1.3) 


— »  /  — ♦  * 

With  this  notation  the  path,  r(a,r^)  ,  is  the  extremum  path  when 


.  .  jL(r»  ,r,a) 
da  5x* 


=  Q 


\  =  x(l,r1) 


xQ  =  *(0,^) 


(A1.4) 


with  similar  equations  in  y(a,r1)  and  z(a,r^)  .  The  fact  that  r(a,r^)  depends 
upon  the  end  points  has  been  explicitly  noted  so  far  as  the  upper  end  point 
is  concerned. 
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Now  consider  the  change  of  9  with  a  change  in  the  value  of  x^  =  x^(l,r^), 
where  a  new  extremum  path  is  chosen  corresponding  to  the  new  value  of  x7  . 


a? 

ax. 


JU 

1  ri  rr 

ah 

X 

dx 

.+  dL 

ax 

j  +  '  do 

aL 

+ 

aL 

ay 

ax» 

axx 

ax 

OX  i 

h 

by1 

axx 

ay 

axx 

1 

r 

da 


ah  az*  ai- 

1z'  bx±  bz 


az 


ax. 


Partially  integrating  and  using  equations  (AlA): 


a? 

ax 

aL 

i 

0=1 

4- 

r 

dy 

aL 

0=1 

4-  - 

axx 

axx 

ax* 

T 

o=0 

ky 

J 

T 

o=0 

at ' 


0=1 


az 


lrr=, 


0=0 


At  the  end  points,  with  a  change  in  only  x^ 


with  all  other 


similar  partials  being  identically  zero.  Thus, 


(A1.5) 
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where  s(r. 

to  (A1.5) 


)  is  the  unit  vector  tangent  to  the  extremum  path  at  the  point  r^  . 
Partial  derivatives  with  respect  to  y^  and  yield  equations  similar 
Consequently ,  the  three  partial  derivatives  combine  in  vector 


notation  to  give: 
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APPENDIX  II.  Evaluation  of  the  Phase  Integral  and  Extremum  Path. 


In  this  appendix,  the  line  integral. 


9  =  I  dr1  *  s(r‘)  f(r* ) 


(A2.l) 


r  (t  ) 

s'  o 


/  1  i  -* 

is  evaluated  to  Ol — r— ]  together  with  the  deviation  ,  §  ,  of  the  extremum  path 

K(D  (  1  \ 

from  the  geometrical  path  to  Ol — g—  !  .  Let: 

\eor  J 


^  =  unit  vector  from  r  (t  )  to  r  , 


s'  o 


§  =  §g  =  vector  difference  between  the  extremum  path  and 


_  A 

geometrical  path,  §  ’  R  =  0  . 


It  is  assumed  that  the  function,  f(r)  ,  is  of  order  unity  and  that 


\vf{r)\  =  . 

\U)  ' 


It  will  he  shown  that 


The  vector  integration  variable  in  equation  (A2.l)  is  given  by: 


r*  =  r‘R  +  §(r’)  ,  0  ^’r*  <:  R  =  jr  -  r  (t  )|  , 

s  o 
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■where,  at  the  end  points  of  the  integration. 


1(0)  =  1(B)  =  0  - 


It  can  he  seen  that 


I^dr*  +  d§ 


|  =  (R  +  lOdr* 


where 


Therefore, 


s(r* ) 


R  +  I1 

fi  +  il'l2 


(A2.2) 


so  that 


dr*  •  s(r‘ )  = 


Vl  +  \%: 


and  equation  (A2.l)  can  he  written  in  the  form. 


=  /dr’  V  1  + 


■*|2  f(r*  ,1(r*)) 


(A2.3) 
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By  definition,  ?  is  an  extremum  if  the  extremum  path,  §(r* ) ,  is  such 


that 


dr 


5fi 


LVi  +  |  §*  i: 


f(r,5) 


f(r,  I) 


k  =  1,2. 


Since,  f  =  0(l)  and  j^fj  =  >  1 5*  |  =  of  ^  )  an&  at)OVe  equation 

'  £0  /  '  <x> 


to  0(l/o>  )  can  he  written  in  the  form: 


dr 


C%(r)  f(r,?)]  =  -g|-  f(r,I)  +  of-ijp)  ,  k  -  1,2 


\  CD  / 


Consequently,  §  should  also  he  of  0(l/o>2)  . 

Assuming  hoth  §  and  §’  are  small,  and  expanding  equation  (A2.k) 
to  0(l/co2)  , 


ir  [§k(r)  f(r>o)]  =  ?) 


+ °  -V) 

?  =  o 


'  CD  1 


(A2.4) 
,  k  =  1,2  . 


2 

Integrating  this  equation  with  respect  to  r  and  re-expanding  to  O(l/co  )  , 


o 


J-  f(r',?(r') 


i  =  o 


(A2.5) 

k  =  1,2. 
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Integrating  this  equation  and  using  the  fact  that  §  at  the  end  points  of  the 
integration  must  he  zero: 


R 

r 


0  = 


r' 

r 


dr‘ 


dr" 


+  R  \(o) 

1  =  0 


or. 


R 

|dr‘ 

t 

i 

/ 

o 


r1 


dr' 


o 


1  =  0  . 


Substituting  this  equation  into  equation  (A2.5), 


r 


o 


R 


I 


|  =  0  o 


-Sr  | 


(A2. 6) 


k  =  1,2 


Finally,  integrating  this  equation, 

(A2.7) 

|  =  0  ,  k  =  1,2, 

2 

which  is  of  0(l/co  ),  justifying  the  above  assumption. 
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To  evaluate  the  phase  integral  S  ,  given  by  equation  (A2.3),  expand 
the  integrand  te  0(l /ot).  Then, 


R 

^  — ♦  — ♦ 

§  =  /  dr*  [f(rI  ,o)  +  — — 2“^ —  f(rJ,o)  +  §  *  Vf(r 


,o)]  +  oL^) 

\  (D  / 


Partially  integrating, 


J 


R 


dr'  |*  •  ?*  f(r',o)  =  -  j  dr1  ?  *  --—--(1^  (r»  ,o)  +  [§  *  i'fCr^o)] 


R 


Therefore,  since  §  =  0  at  the  end  points. 


R 


§  = 


[f(r’,o)  +  I 


(vf(r\o)  - 


d 

dr* 


(f!*))]  +  0 


o 


Partially  integrating  once  more  after  substituting  from  equations  (A2.4),  the 
final  expression  becomes: 


4 


* 


(A2.8) 


where  the 


V  (r’)  are  given  to  sufficient  order  by  equation  (A2 .6). 
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